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Magnetism and anti-ferroelectricity in MgBg 
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We report on a density functional theory study demonstrating the coexistence of weak ferromag- 
netism and anti-ferroelectricity in boron-deficient MgBg. A boron vacancy produces an almost one 
dimensional extended molecular orbital, which is responsible for the magnetic moment formation. 
Then, long range magnetic order can emerge from the overlap of such orbitals above percolation 
threshold. Although there is a finite density of states at the Fermi level, the localized nature of 
the charge density causes an inefficient electron screening. We find that the Mg^"^ ions can displace 
from the center of their cubic cage, thus generating electrical dipoles. In the ground state these 
order in an anti-ferroelectric configuration. If proved experimentally, this will be the first material 
without d or f electrons displaying the coexistence of magnetic and electric order. 

PACS numbers: 81.05.Zx, 71.55.-i, 75.10.Lp, 75.30.Hx, 77.80.-e 
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The discovery of weak ferromagnetism in La-doped 
CaBg by Young et al. [l[ ignited a surge of research on 
alkaline-earth hexaborides, XBg (X — Ca, Sr, Ba) [2|-[8[. 
Initially the origin of the magnetism in these materials 
was attributed to possible transition- metal impurities [9| . 
However, subsequent experiments in high-purity CaBg, 
SrBg and BaBg crystals 10l4l2J clearly indicated that the 
magnetism may be an intrinsic property of the XBg fam- 
ily, and not due to contaminations. Theoretical predic- 
tions for defective XBg assigned to boron vacancies, Vb, 
the primary cause of the magnetism. These introduce flat 
bands near the Fermi level, which spin-split due to the 
strong exchange interaction in the B p-shell [3|, |j, [iSl, \\A \ . 
Such predictions have been recently supported by exper- 
imental evidence [sl, lla |l6[ , which also extends to CaBg 
and SrBg thin films [7| . 

In contrast to the heavy compounds of the XBg family, 
MgBg has received only little attention. This is mainly 
due to an early theoretical calculation, performed with 
empirical methods, suggesting the poor stability of MgBg 
compared to the more thermodynamically stable MgB2, 
MgB4, MgBr and MgBi2 [i3|- In addition the small Mg 
ionic radius, second only to Be among the alkalines, may 
result in a high ionic mobility and thus in Mg segrega- 
tion. However, in contrast to this early prediction MgBg 
has been often observed experimentally [18[, most no- 
tably with conclusive high resolution tunnelling electron 
microscopy [l9|. In this letter we report on a density 
functional theory (DFT) study of the structural, elec- 
tronic and magnetic properties of Ve-rich MgBg. Our 
results point to a ferromagnetic ground state, which co- 
exists with an anti-ferroelectric dipolar order. As such 
MgBg appears to be the first example displaying the co- 
existence of magnetic and electric order in a material not 
containing ions with partially occupied d or / shells. 

Calculations are performed with the DFT numerical 
implementation contained in the Siesta code [20|. The 
Perdew-Burke-Ernzerhof form of the generalized gradi- 



ent approximation |2l| is used for the structural op- 
timization and for the calculation of the phonon dis- 
persion. Additional calculations of the electronic prop- 
erties use the Perdew-Zunger local density approxima- 
tion (LDA) [2a|including atomic self-interaction correc- 
tions (ASIC) [23|. The core electrons are described by 
norm-conserving pseudopotentials with partial core cor- 
rections j24j , while the wave- function is expanded over a 
double-zeta numerical atomic orbital basis set including 
the following orbitals: Mg 3s, Mg 3p, B 2s and B 2p. The 
charge density and the DFT potential are determined 
over a real-space grid with an equivalent mesh cutoff en- 
ergy of 350 Ry, but a value of 700 Ry is used for calcu- 
lating the phonon spectra. A 1728 k-points dense mesh 
samples the Brillouin zone. Atomic relaxation is per- 
formed by conjugate gradient until the forces are smaller 
than 0.04 eV/A (0.001 eV/A for the phonon spectrum). 
Phonons are calculated by using the elementary unit cell, 
hence only optical modes are considered. 

We start our analysis by first looking at the defect-free 
crystal. XBg (X=Be, Mg, Ca) crystallizes in a halite cu- 
bic structure with the two inequivalent sites occupied re- 
spectively by the divalent X^+ ion and a Bg~ octahedral 
cluster [see Fig. [T{a)]. The lattice parameter is deter- 
mined by the size of the Bg clusters and it is calculated 
to be 4.115 A for MgBg. The position of the cation how- 
ever depends on the specific element. After structural 
relaxation Ca remains in the center-symmetric position, 
while Mg moves towards the face of the cubic cage and 
nests at 1.143 A away from it. Interestingly, Be moves 
completely into the Bg plane. 

The potential energy surface of Mg^"*" is plotted in 
Fig. IHb) over one of the median planes cutting across 
the cubic cage, with the zero energy taken at the center- 
symmetric position. We observe four equivalent minima 
in the plane, which translate into six minima over the cu- 
bic volume. The potential barriers between the minima 
are around 70 meV, which indicate that the displaced 
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FIG. 2: (colour online) Electronic structure of defect-free 
MgBe. In (a) we present the type-C antiferroelectric order 
for a 3 X 3 X 3 supercell. The arrows indicate the direction 
of the electric dipoles (atoms and bonds are removed for sim- 
plicity). In (b) a comparison of the electronic band structures 
for the center-symmetric and the type-C anti-ferroelectric or- 
dered system obtained for 2x2x2 supercell. The correspond- 
ing PDOS is presented in the rightmost panel. 



FIG. 1: (colour online) MgBg structural properties. In (a) 
the XBe (X=Be, Mg, Ca) crystal structure, where the gray 
volume includes the unit cell. Note the different displacement 
of Ca, Mg and Be within the Be-defined cubic cage, (b) Po- 
tential energy surface (isovalues are in units of meV) for Mg^+ 
across the median plane of the cubic cage (shown in the in- 
set). Phonon band structure for the center-symmetric MgBe 
crystal (c) and center-symmetric CaBe crystal (d). 



Mg positions arc stable at room temperature. Addi- 
tional evidence of the Mg displacement is provided by 
the MgBg phonon spectrum [Fig. [Ijc)] for the center- 
symmetric structure, which presents a triply-degenerate 
negative frequency mode. Note that the same is absent 
for the CaBg crystal [Fig. [Ijd)], which remains undis- 
torted. 

Given the ionic character of MgBg, the displacement 
of Mg^+ generates an electrical dipole parallel to the 
direction of the displacement. Standard (LDA) Berry- 
phase calculations give us an electrical polarization of 
57.6 /^C/cm^ for an hypothetical ferroelectric phase. This 
however does not represent the correct long-range ferroic 
order. The overall dipolar arrangement is then computed 
by comparing total energy calculations for a 2 x 2 x 2 su- 
percell with different dipole configurations (ferroelectric, 
antiferroelectric of type A and type C and ferroelectric 
with only one dipole chain rotated) l25[ . We find the en- 
ergy minimum at the type-C anti-ferroelectric arrange- 
ment, which is illustrated in Fig. ^a.). Such a state is 
226 meV lower than the ferroelectric one (per elementary 
unit cell), which suggests that it might be a relatively 
high-temperature phase. 

The projected density of states (PDOS) of MgBg in- 
dicates that the material is semiconducting with a rel- 
atively small indirect (LDA) band gap of about 0.5 eV 
[see Fig. [H^b)]. This reduces to 0.27 eV if the calcula- 
tion is performed at the GGA level. The conduction and 



valence bands have respectively Mg-s and B-p character, 
as expected from its ionic Mg^+(B6)^~ nature. An ionic 
bonding favours the Mg^"'' ion to move off-center in order 
to minimize the electrostatic energy. The same mecha- 
nism is ineffective for the heavier CaBg simply because 
of the much larger ionic radius of Ca^+. Furthermore, 
as there is little covalent bonding between Mg and Bg 
also the interaction between the dipoles has electrostatic 
nature and so an anti-ferroelectric order is expected [26| . 
Finally we note that the Mg displacement has little ef- 
fect on the band structure as the anti-ferroelectric ground 
state remains insulating with an LDA gap of 0.73 eV 
[Fig.Hb)]. 

Having characterized the electronic structure of defect- 
free MgBg, we now move to discuss the effect of Vb in the 
crystal. The MgBg crystal bares symmorphic symmetry, 
so that its bands can be directly related to the energy lev- 
els of the Bg cluster. The molecular states [Fig. ^h)] are 
characterized according to the irreducible representations 
of the O/i group. The low lying Tin state is of special 
interest. In fact a vacancy in Bg crystal moves the Ti„- 
derived level upwards into the pristine gap of the crystal. 
Ti„ is threefold degenerate with each of the levels in the 
manifold dominated by only one of B-p^ {i ~ x, y, z) or- 
bitals. A vacancy also removes an electron so that B5 is 
left with a single electron in non-degenerate orbital level. 
We then expect magnetic moment formation [27[ with 
the spin-polarized charge density distributing along the 
direction where the B atom has been removed. 

This qualitative picture is further confirmed with nu- 
merical DFT calculations. Since the self-interaction er- 
ror in LDA usually precludes the correct description of 
defects in semiconductors, here we employ the ASIC 
scheme [23||, which has been already proved successful 
in describing anion defects in the closely related MgO in- 
sulator [28| . In particular we perform calculations by us- 
ing different values of the screening parameter a, ranging 
from a={) (LDA) to a=l (full ASIC), and compare spin- 
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FIG. 4: (colour online) Local density of states around the 
Fermi level of a MgBg crystal with one (a) and two (c) vacan- 
cies per 2x2x2 supercell. In the panels (b) and (d) we show 
the DOS projected over different sub-set of B atoms. In par- 
ticular we project over B atoms along the directions indicated 
by the arrows in the panels (a) and (c). 



FIG. 3: (colour online) Electronic structure of Vs-rich MgBg. 
In (a) the geometry of a supercell containing one Vb. (b) 
shows the energy state diagram of a free-standing Bg cluster. 
The red arrow indicates the re-arrangement of states upon 
introduction of a vacancy in MgBe crystal. The electronic 
band-structure and DOS for the cell containing one vacancy 
are in panel (c). 



polarized (SP) and non-spin-polarized (NSP) results. 

Let us first discuss the magnetic moment formation 
associated to a single Vb in one 2x2x2 supercell [see 
Fig.El^a)]. LDA yields a cell magnetic moment of 3 m/ie. 
Importantly such a moment increases in the ASIC cal- 
culations. In fact we find 414 ni/^B, 509 TafiB, 565 m^ue 
and 617 m/.tB respectively for a = 0.25, 0.5, 0.75 and 
1. Also the total energy difference between the SP and 
the NSP solutions, AE{a) = Esp — -Bnsp, varies with 
the degree of ASIC considered (a). In particular we find 
AE= 2 meV for LDA and -23 meV, -44 mcV, -116 mcV 
and -296 meV respectively for ASIC with a = 0.25, 0.5, 
0.75 and 1. As such we find that both the magnetic mo- 
ment and its energy stability increase with increasing the 
amount of self-interaction corrections included. A value 
a = 0.5 is usually considered appropriate for semiconduc- 
tors so that we estimate a magnetic moment of around 
1/2 /iB, with a formation energy well above the room 
temperature. 

A deeper insight into the origin of the magnetic mo- 
ment is obtained by looking at the band structure of the 
supercell containing one Vb [Fig.|31Jc)]. As expected from 
the discussion about Bg now the system has a metallic 



electronic structure with three bands crossing the Fermi 
level, ^F [labelled with 1, 2, and 3 in Fig. [Sljc)]. Their 
dispersion is typical of a semi-metal with both electron 
(from M to X and at R) and hole (at F) pockets at the Ep. 
As the three bands are non-orbital-degenerate the system 
is not prone to the Jahn- Teller distortion, which might 
compete with the formation of the moment. Importantly 
only the electron carrying band [3 in Fig. [3Kc)] spin- 
splits, while the other two remain non-spin-polarized. 
Recalling the fact that the Stoner exchange parameter 
is rather large in the 2p shell 
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the spin-splitting oc- 
curs in bands with a narrow band-width, i.e. in bands 
that can provide a large DOS at Ep. This is the case of 
the electron carrying band (3), which has a band-width 
of 0.68 cV to be compared with 1.3 eV and 1.2 eV for the 
other two. Such a picture of magnetic moment formation 
due to narrow bands is robust with both the number and 
the position of the vacancies as we tested for a 2 x 2 x 2 
supercell with two Vb's at different positions. 

The unique directional nature of the VB-derivcd elec- 
tronic bands at Ep is provided in Fig. |4l[a), where we 
show the local density of states (LDOS), obtained within 
the energy range (£'F-0.05eV, i?F-|-0.05cV). From the fig- 
ure it emerges that the electron cloud is delocalized only 
along the direction defined by the center of vacant B5 
cluster and by the position of Vb , as expected from the 
Ti„ symmetry. The contribution to the DOS from the B 
atoms which do not lie along such a preferential direction 
is negligible [Fig. IUb)]. The presence of multiple Vb's 
does not alter this picture. Fig. UJ^c) shows the LDOS 



of two Vb's at different positions, with formation of two 
directional states oriented in different directions and not 
overlapping among each other. 

The formation of a moment is not a sufficient condi- 



tion for long range magnetism |28| . The picture emerging 
from our electronic structure calculations is that of mag- 
netic moments associated to directional extended molec- 
ular orbitals at the B5 clusters. We then expect mag- 
netism only when the magnetic-moment-carrying levels 
of different B5 clusters overlap. We have checked that 
distant Vb 's are magnetically coupled by calculating the 
total energy difference between a ferromagnetic and an 
antiferromagnetic configuration of a 4 x 2 x 2 supercell. 
We find that, when the two B5 clusters are 8.23 A apart, 
there is an energy difference of 32 meV/vacancy (a = 0.5) 
in favor of the ferromagnetic solution. Although more 
calculations are needed at other Vb concentrations and 
for different Vb-Vb separations, we may yet conclude 
that our calculations point to both the moment forma- 
tion and exchange coupling between distant moments. A 
long range order is then expected when, typically ran- 
domly distributed, the magnetic moments form a perco- 
lation network. Two important considerations must be 
made here. First, as one missing B over the Bg cluster 
is responsible for the moment and percolation is between 
the Bg clusters, a cluster percolation threshold of x is 
achieved by removing only x/6 B ions. Second, we re- 
mark that the magnetism is related to the Bg sub-lattice 
only, so that our results should be transferable to other 
XBg materials. However, we expect that the local electri- 
cal dipole associated to the Mg^+ ions in cages bordered 
by B5 clusters, will be different in direction and magni- 
tude from that of the vacancy-free case. 

Finally we wish to spend a few words on the possi- 
bility of the coexistence of both the ferromagnetic and 
anti-ferroelectric order. In general a ferroic state cannot 
be sustained in a metal because of the screening from the 
free electrons. The situation here is however more com- 
plex as, although there is density of states at the Fermi 
level, this is associated to partially localized states. Their 
screening can be inefficient to a point that local electric 
dipoles still form [29|,|30[. Our results indeed confirm such 
a hypothesis as the Mg ions remain in their off-center po- 
sitions as B vacancies are introduced in the cell. Notably, 
the strong directional nature of the charge density at the 
Fermi level suggests that such an inefficient screening per- 
sists locally even above the percolation threshold. 

In conclusion, we have investigated structural, elec- 
tronic and magnetic properties of ideal and defective 
MgBg-s crystals. These present the coexistence of weak 
magnetism and potentially anti-ferroelectricity. The lo- 
cal magnetic moments originate from a VB-related impu- 
rity band and the dipole moments from the displacement 
of the Mg^+ ions within the cubic cage. Long-range mag- 
netism can be sustained only above the magnetic percola- 
tion threshold and the anti-ferroelectric state can coexist 



with a finite DOS at the Fermi level because of the in- 
complete screening from the localized charge. If proved 
experimentally, this will be the first example of coexis- 
tence of magnetic and electric order in a material not 
incorporating ions with partially occupied d or f shells. 
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